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ABSTRACT 

Sets  of  systems  of  ordinary  functional-differential  equations  with 
Volterra  type  functionals  under  sets  of  initial  values  are  considered.  Upper 
and  lower  bounds  are  constructed  for  the  sets  of  all  solutions.  Classes 
of  such  problems  are  given  where  these  bounds  are  optimal.  The  main  tool 
is  a Lemma  of  Max  Muller  on  inequalities.  Also  ideas  from  interval 
mathematics  are  used. 

AMS  (MOS)  Subject  Classification:  34J05 

Key  Words:  Functional-differential  equations  and  inequalities , Volterra  type. 
Solution  set.  Optimal  bounds.  Maximal  interval  solution 

Work  Unit  Number  1 (Applied  Analysis) 

SIGNIFICANCE  AND  EXPLANATION 

If  differential  equations 

u'(t)  = f(t,u(t)),  u (0)  = a 

appear  in  Applied  Mathematics  there  is  normally  not  just  one  right  hand  side 
f(t,u).  Instead  of  this  a whole  set  {f}  of  right  hand  sides  must  be 
considered.  This  is  due  to  many  facts  such  as:  data  errors,  data  intervals 
obtained  from  measurements,  approximation  of  f by  a more  suitable  function, 
poor  knowledge  of  the  laws  involved  etc.  The  same  is  true  for  the  initial 
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"value"  ct  which  is  usually  a set  {a}.  Hence  the  above  initial  value 
problem  has  to  be  replaced  by  the  inclusion  problem 

u' (t)  e {f (t,u(t))},  u (0)  e {a}  . 

It  is  normally  completely  impossible  to  solve  all  the  real  problems 
which  are  combined  in  this  set  of  problems.  The  goal  of  the  following  paper 
is  therefore  to  find  at  least  lower  and  upper  bounds  to  the  set  of  all  such 
solutions.  This  can  always  be  done.  Since  these  bounds  are  sometimes  very 
pessimistic,  classes  of  such  problems  are  given  where  the  bounds  obtained 
are  optimal . 

The  main  ideas  of  this  paper  are  also  valid  in  the  more  general  case 
where  f does  depend  as  a functional  upon  the  unknown  solution  u.  This 
is  written  in  the  form  f (t ,u (t) ,u ( • ) ) . Therefore  the  theory  of  this  paper 
also  includes  integro-dif ferential  equations  and  difference-differential 
equations.  Sets  of  such  problems  do  occur  for  example  in  Economics  and  in 
Biology. 


BOUNDS  FOR  THE  SET  OF  SOLUTIONS  OF  FUNCTIONAL-DIFFERENTIAL  EQUATIONS 


Karl  L.  Nickel* 


1 . Introduction 

In  the  following  paper  systems  of  functional-differential  equations 

(1)  u'(t)  = f (t.u(t)  ,u(0)  for  0 < t < T 

s 

are  considered  under  the  initial  corelitions 

(2)  u(0)  = a . 

Herein  u = (u  ,u  , — ,u  ),  f = (f  ,f  ,...,f  ) and  a = ( a , ,a_ , . . . ,a  ) are  n-vectors. 
12  n 12  n 12  n 

As  usual  u or  u(t)  means  the  value  of  the  function  u at  the  point  t;  moreover 
u‘ (t)  = du/dt.  Opposite  to  this  the  notation  u(-)  means  that  u is  regarded  as  an 
element  of  the  class  of  admissible  functions.  Hence  f(-,-,u(-))  is  a functional  on  u; 
in  what  follows  only  special  "Volterra"  functionals  will  be  regarded. 

If  f is  continuous  then  the  system  (1),  (2)  is  equivalent  to  the  system  of  func- 
tional-integral equations 

t 

(3)  u(t)  = a + / f (s,u (s) ,u ( • ) )ds  for  0 < t < T . 

0 

It  is  the  subject  of  the  following  paper  to  find  bounding  functions  v(t),  w(t)  such 
that  for  every  solution  u(t)  of  (1),  (2)  or  (3) 

(4)  v(t)  < u(t)  < w(t)  for  0 < t < T . 

s & s s 

Hence  the  classical  theory  of  maximal  and  minimal  solutions  for  differential  equations 
appears  as  a special  case  of  these  results. 

If  a solution  u of  (1),  (2)  is  uniquely  determined  then  it  is  trivial  that  (4)  is 
satisfied  for  w :=  v :=  u.  It  is  therefore  interesting  to  switch  to  a more  general 
problem:  Let  {a}  be  a set  of  initial  values  and  let  {f}  be  a set  of  right  hand  sides 
to  (1) . Then  the  more  general  initial  value  problem 
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(5) 


u'  { t)  e {f  (t,u(t)  ,u(-)) } for  0 < t < T , 

(6)  u (0)  * {a} 

is  considered.  Let  {u}  = {u(t)}  be  the  set  of  all  solutions  of  (5),  (6)  Again  two 
functions  v(t),  w(t)  are  looked  for  such  that  (4)  is  true  for  any  solution  u e {u}. 

If  one  writes  [v,w]  for  the  function  interval  from  the  two  bound  functions  v and  w 
then  this  can  be  written  as 

(7)  {u}  t [v,w]  . 

It  is  in  general  quite  simple  to  find  rough  bounds  v,w.  In  what  follows  special 
emphasis  is  therefore  given  to  the  look  for  "optimal"  bounds.  Here  "optimality"  means 

the  following:  let  there  exist  the  infimum  and  the  supremum  of  the  set  {u}  such  that 

» 

(8)  v = inf{u),  w = sup{u}  . 

In  that  case  one  can  call  (v,w)  the  "interval  hull"  of  {u).  It  is  the  goal  of  this 
paper  to  find  classes  of  sets  { f } and  {a}  such  that  (7)  and  (8)  are  true. 

In  order  to  get  such  results  a lemma  of  Max  Muller  (1927)  on  differential  inequalities 
is  essential.  This  lemma  has  been  published  exactly  50  years  ago.  For  decades  however, 
it  remained  widely  unnoticed.  In  what  follows  this  lemma  will  be  extended  to  the  case 
of  functional-differential  inequalities.  This  will  be  done  by  extending  an  old  paper 
of  the  author  (Nickel  (1961)). 

It  should  finally  be  remarked  that  the  problem  of  this  paper  and  some  of  the  formula- 
tions have  been  strongly  influenced  by  the  ideas  of  interval  mathematics. 
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2 . Notations  and  assumptions 


Let  n t W,  0 < T f K,  I :=  [0,T]  , IQ  :=  (0,TJ,  a e Rn.  The  n-vectors  a,  u(t), 

f(t,u,u(*))  are  written  as  a = (a,,a_,...,a  ),  u = (u,,u. ,u  ),  f = (f ,, f f ). 

12  n 12  n 12  n 

Together  with  the  k-th  component  of  the  vector  u also  the  n-l-vector 

u = (u,  ,u,,...,u  ,>u.  u ) is  used, 

k 12  k-1  k+1  n 

Let  the  class  2 of  the  admissible  functions  be  defined  as  the  set  of  all  function 
vectors  u : I -*  Kn,  continuous  on  I such  that  the  derivative  u'  = du/dt  exists  in 

The  notation  u(-)  means  that  u is  an  element  of  the  class  2.  Opposite  to  this 
u,  u'  or  u(t),  u'(t)  mean  the  values  of  these  functions  at  the  point  t. 

Inequalities  v(t)  < w(t)  or  v < w are  always  meant  componentwise  as  v^(t)  1 wk<t) 
for  k = 1(1)  n.  Inequalities  of  the  kind  v(-)  < w(-)  are  meant  both  componentwise  and 
pointwise  for  all  points  in  the  definition  set. 

For  v < w the  interval  [v,w]  is  defined  as  the  set  (v,w]  = [v(t),w(t)l  :=  {z  f k'1 
v(t)  < z < w(t)>.  Similarly  lv(*).w(-)l  :=  {z  e Z | v ( - ) < z ( - ) < w ( - ) } . 

Let  the  dependence  of  f^(t,u,u(*))  of  any  component  ^or  k = l(l)n  be 

that  of  a "Volterra"  functional.  Here  a functional  g(t,z(-))  is  called  "Volterra" 
if  it  represents  a mapping  in  ]R  such  that  the  value  of  g at  the  point  t depends 
only  on  the  value  z(s)  for  0 < s < t (g  depends  only  upon  “the  past”  of  the  function 
z(-)). 

Examples  of  Volterra  functionals  are 


t 

9 (t,2  ( • ) ) 

:=  / K(t,s,z(s))ds  , 

0 

(9) 

g(t,z(-)) 

:=  z (t  « t)  with 

0 < T 

r—i 

V II 

C z(t  - s)  for 

0 < s 

< t , 

(10) 

g(t,z(*)) 

= = ( 

= 

e 

a (t  - s)  for 

o 

II  A 

ft 

< S 

with  some  given  : 

function 

a(t)  for 

-s  < t < 0 . 

These  examples  show  that  the  theory  given  in  this  paper  can  be  applied  to:  differential 
equations,  (Volterra)  integro-dif ferential  equations,  difference-differential  equations 
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with  retarded  argument  and  naturally  also  to  combinations  of  these  equations. 

In  order  to  simplify  the  results  the  following  notation  wi‘11  be  used  for  the 

components  f^  of  f(t,u,u(*)):  the  second  argument  u f ]Rn  of  f is  broken  up  in 

the  component  u with  the  same  index  k as  f,  and  in  the  rest  vector  .u: 

k k k 

fk  = fk*t,uk'ku’u*'"  * 

Furthermore  it  is  suitable  to  have  a special  notation  for  the  set  of  all  functions  f^_ 
if  the  arguments  lie  in  certain  intervals.  This  will  be  denoted  by 

'fk<t,Uk'  *kV'kWl  ' ) ) :=  (f  (t,uk,ku,u(-))  lku  e Ij^v.^w]  ,u  ( • ) e [v(-)  ,w(-)J  } . 

From  section  7 to  the  rest  of  this  paper  only  functions  f will  be  regarded  which  are 
partially  monotone.  The  corresponding  definitions  will  be  given  in  section  7. 


-4- 


3.  Existence 


— 

The  following  theorem  is  the  extension  of  the  well  known  Peano  existence  theorem 
for  systems  of  differential  equations: 

Theorem : 

I 

f be  defined  and  continuous  on  I x Kn  x Z.  Then  (1),  (2)  and  (3)  are 
equivalent.  If_  f is  bounded  there  exists  (at  least)  one  solution  ue  Z of  (1),  (2). 
— f is  not  bounded  then  there  exists  a solution  of  (1),  (2)  at  least  in  a largest 
interval  0 < t < t < T. 

Proof . The  equivalence  is  trivial.  For  the  existence  the  fixed  point  theorem  of  Schauder 
is  applied  to  equation  (3) . The  main  ideas  are  exactly  the  same  as  in  the  case  of 
differential  equations.  They  are  described  in  the  book  of  Walter  (1970),  p.  23-25. 

j 

j' 


j 


4 . The  Lemma  of  Max  Muller 
Lemma : 


Let  the  functions  v,w  f Z with  v < w satisfy  the  following  inequalities 


(11) 

v (0) 

< a < w(0)  , 

(12) 

vk  < {fk(t,V 

[kV'kW'  ' lv<‘) 

(13) 

wk  ” {fk(t'v' 

[kV,  W]  , [v  ( •)  ,w(  •)  ] } 

iSSL  * e J0 

and  k = 1 (1)  n . 

Then  any  solution  u e Z o£  (1)  , (2)  is  bounded  by 

(14)  v(t)  < u(t)  < w(t)  for  tel. 


Corollaries : 


1 )  If  u e Z is  a solution  of  the  inequalities  u(0)  £ a,  u'  < f ( t , u , u ( • ) ) in 


1q  then 


u(t)  < w(t)  for  tel. 


2)  Similar 

u (t)  > v (t)  for  tel 

for  any  solution  u e Z of  the  inequalities  u(0)  > a,  u'  > f(t,u,u(*))  ill  1^. 

3)  If  all  f are  strictly  monotone  (increasing  or  decreasing)  with  respect  to 
(at  least)  one  of  the  components  of  u("),  it  then  suffices  to  have  the  > - and 
<-signs  in  (12)  and  (13)  instead  of  the  >-  and  <-signs . 

Remarks : 

1)  This  Lemma  has  been  formulated  and  proven  by  Max  Muller  (1927)  as  Theorem  5 
on  the  pages  13  to  15  for  the  special  case  where  f does  not  depend  upon  u(‘).  See 
W.  Walter  (1970),  p.  93-94. 

2)  The  original  notation  of  M.  Muller  was  very  inconvenient.  It  has  here  been 
replaced  by  the  interval  notation. 

3)  Kindly  note  that  there  are  no  assumptions  to  be  made  with  respect  to  the  function 
f,  such  as  continuity,  monotonicity  etc. 

Proof : The  original  proof  of  M.  Muller  carries  right  over  to  this  case. 


j 
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5.  Example 


In.  the  following  example  n = 2.  The  functional  (9)  is  used  for  t = 1/2  together 
with  the  other  functional 

t , , 

/ (u  (s)  + u (s)  Ids  . 

0 1 * 

Let 

fjlt.u.uf.))  : = -2u2  + u2/(l  + t2)  + (ui(t/2)  + u2<t/2))/2  , 

f2(t,u,u(-))  :**  1 - sin(Tiu^/2)  + 2u2<1  - 2u2>  + / (u2(s)  + u2(s))ds  , 


° < °1'  °2  < 1 ‘ 

Define  = v 2 :=  0,  w^  = w2  :=  1 for  t > 0.  Then  one  verifies  easily  for  t > 0 
fjlt.Vj. Iv2,W2],lvi(-).W1(-)].tV2(-).W2(-)])  = 0 4 [0,1]/(1  4 t2)  4 [0,1]  > 0 , 
f1(t,w1,[v2,w2l,|v  (.),W  (•)],[v2(-),w2(-)))  =-24  [0,1]/(1  4 t2)  + (0,1)  < 0 , 
f2(t,V2,...)  = [0,1]  404  [0,1]  > 0 , 
f2(t,w2,...)  = [0,1]  - 2 4 [0,1]  < 0 . 

Since  f^  and  f2  are  strictly  monotone  increasing  with  respect  to  both  components 
u^f-)  and  u2 ( - ) the  third  corollary  to  the  lemma  can  be  used.  This  gives  the  a priori 
estimate  0 < u^(t),  u2(t)  < 1 for  t > 0 for  any  solution  u. 
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6.  Uniqueness  and  error  bounds 

For  ordinary  differential  equations  all  uniqueness  theorems  can  be  derived  from 
the  Lemma  of  M.  Muller.  In  his  original  paper  he  did  use  the  lemma  exactly  for  that 
purpose.  Similarly  probably  all  known  a posteriori  bounds  can  be  proven  with  its  help. 
Therefore  a very  large  number  of  such  theorems  can  immediately  be  proven  also  for  the 
case  of  functional-differential  equations.  It  is  possible  to  translate  all  the  results 
from  the  book  of  W.  Walter  to  these  extended  systems  of  equations.  Some  first  results 
have  already  been  given  in  the  old  paper  by  K.  Nickel  (1961) . It  is  however,  not  the 
purpose  of  this  paper  to  publish  such  theorems. 


7 . Monotonicity  conditions 

Definition  (unconditional!/  partially  isotone/antitone/monotone) : 

Let  g (Xj  ,Xj » . . . , JO  be  a mapping  g : D -*•  R with  D C Rm.  Again  the  notation 
g(x^,...,x^)  = 9(x^.^x)  is  used.  The  function  g is  called  unconditionally  partially 
iso tone  or  antitone  on  D with  respect  to  the  variable  x^  if 

g(Vkx)  ^ 9(Vkx)  for  yk  = xk  or  yk  i \ 

and  for  all  Cx^.^x)  , (yk<kx)  * D.  A function  which  is  either  unconditionally  partially 
isotone  or  unc.  part,  antitone  is  called  unconditionally  partially  monotone. 

Kindly  note  that  the  function  gfx^.x^)  :=  x^-x2  *s  “"conditionally  partially 

2 

monotone  with  respect  to  and  x2  on  D :=  [0,“)  x (0,“),  but  not  on  D R . 

Definition  (monotonicity  class  M)  : 

let  the  class  M consist  on  all  functions  f(t,u,v(-))  for  which  the  following 
is  true:  Each  function  f (t .u^.^u.v  ( • ) ) is  unconditionally  partially  monotone  on 
I0  * * - 11°  x Z with  respect  to  any  component  of  ^u  and  v(-),  but  not  necessarily 
with  respect  to  u^. 

If  f e M then  it  is  convenient  to  write 

fk(t,uk,ku,u< •) ) = fk(t,uk,kut,kul,u(*)+,u(-)l)  . 

This  clearly  means  that  the  vectors  ku  and  u(-)  are  divided  in  the  two  sets  of 
components  for  which  f is  isotone  (t)  and  antitone  (!) . 

Let  ft  M.  Then  in  the  Lemma  of  M.  Muller  the  (rather  inconvenient)  inequalities 
of  sets  can  be  replaced  by  real  inequalities  (which  are  much  simpler  to  handle)  . In 
this  case  the  Lemma  reduces  to  the 
Special  case  of  the  Lemma: 

Let  f t M.  Assume  that  (11)  is  true  for  functions  v,w  e Z with  v < w and 
that  furthermore 

vk  < fk(t,vk,kv+,kwi,v(-)+,w(-)+)  , 
w£  > fk(t,wk,kwt,kv+,w(-)+,v(.)+) 

for  t t 1^  and  k = l(l)n  . 
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Then  any  solution  u t Z of  (1).  (2)  is  bounded  by  (14);  furthermore  the  corollaries 
1)  to  3)  hold. 

The  proof  comes  immediately  from  the  Lemma  together  with  the  definition  of  the 
class  M. 


8.  Construction  of  bounds 

In  what  follows  the  special  case  of  the  Lemma  will  be  used  to  construct  bound  func- 
tions v and  w.  The  above  problem  (1) , (2) , (4)  will  be  replaced  by  a more  general 
since  no  additional  difficulties  are  generated  by  that  extension: 

The  initial  data  a in  the  initial  conditions  (2)  are  to  be  replaced  by  a set  la}. 
For  simplicity  assume  {a}  C [a, a]  with  two  vectors  a, a f Rn  and  a < a such  that 
(15)  a, a e {a}  C [a,a]  . 

Hence  the  initial  conditions  (1)  are  to  be  replaced  by  (6).  Let  lu)  be  the  set  of 
all  solutions  of  (1)  under  the  set  of  all  initial  conditions  (6) . Wanted  are  bounds 
v,w  such  that  (7)  is  true. 

Assume  that  f is  continuous  and  bounded  on  I x Rn  x Z . The  following  sequence 
of  problems  for  v « R is  considered: 


(16) 


Vk  = fk(t'Vk'kVt'kWl'kV(')t'kW(-)*)  ' v ' 

Wk  = fk(t'Wk'kW+'kW'kW(-)t'kV(,)l)  + i 
for  t e IQ  and  k = l(l)n  , 

(17)  v(0)  = a - 1/v,  w(0)  = a + 1/v  . 

By  the  existence  theorem  of  section  3 there  exists  for  every  v e R at  least  one 

solution  v,w  t Z of  the  coupled  system  (16),  (17).  An  arbitrary  solution  is  picked 

V V 

up  and  called  (v  ,w  ) . Then 

v v+l  - v+1  v 

v < « < u < w <w 

for  all  v e R and  for  any  solution  u t Z of  (1),  (2).  This  follows  by  the  special 

case  of  the  Lemma  and  by  the  definition  of  the  right  hand  sides  in  (16)  and  in  (17). 

ice  the  sequences  {vV}  and  (wV)  are  monotone  and  bounded.  As  sequences  of 

measurable  functions  they  have  measurable  limit  functions  v(t)  :=  sup{vV(t)}  and 
— V 

w(t)  :=  inflw  (t)}.  Furthermore 

(18)  v < u < w . 

One  can  show  as  usual  that  these  sequences  are  uniformly  convergent  on  I and  that 
v,w  f Z (see  W.  Walter  (1970),  p.  68).  Furthermore  the  pair  (v,w)  satisfies  the 
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following  functional-differential  system  in  I0  consisting  of  2n  equations 


(19) 


*»k  “ fk(t.wk,kw+,kvi,w(.)t,vC)4) 


under  the  2n  initial  conditions 


(20)  v(0)  = a,  w(0)  = a . 

If  u e Z is  a solution  of  (1),  (2),  then  the  pair  (u,u)  is  also  a solution  of  (19), 
(2).  In  general  the  reverse  is  however,  not  true,  i.e.  the  functions  v,  w are  in 
general  no  solutions  of  (1) , (6) . Hence  they  are  not  minimal  or  maximal  solutions  to 
(1)  , (6)  in  the  usual  sense  (see  however,  section  11).  In  what  follows  the  interval 
Iv.w)  will  be  called  maximal  interval  solution  of  (1),  (6). 

The  reason  for  this  notation  comes  from  the  following:  Define  the  interval  operator 
F|v,w)  = (F^ . . . ,Fn>  by  its  k-th  component  as  follows: 


Fk(v,w)(t)  :=  (ak  + / fk(s,vk,kv+,kw+,v(')+,w(-)+)ds  , 


“k  + / fk<s,wk,kw+,kv+,w(-)  + ,v(-)+)ds) 


Any  solution  u of  (1) , (6)  satisfies 

u e F lu ,u]  . 

Moreover  the  interval  (v,wj  is  a fixed  interval  of  the  operator  F by  (19)  and  (20) . 

By  construction  [v,w]  is  the  smallest  fixed  interval  of  F for  which  (18)  is  true. 

If  a = a = a and  if  maximal  and  minimal  solutions  of  (1) , (2)  exist  then  they  are 
equal  to  v and  w. 

This  idea  consists  therefore  in  replacing  the  usual  ordering  relation  < (component- 
wise with  respect  to  k and  pointwise  with  respect  to  t)  by  the  inclusion  C as  a 
new  ordering  relation  (also  componentwise  and  pointwise) . 
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9 . Bounds  for  the  solutions  of  sets  of  functional-differential  equations 


The  equations  (1),  (2)  are  now  being  replaced  by  inclusions  (5),  (6). 

Theorem : 

Let  (15)  be  satisfied  by  the  set  of  initial  values  {a}.  Assume  that  there  exist 
two  right  hand  sides  ft  (f)  with  f < f such  that  for  all  solutions  u t Z of_  (5)  , 

f (t,u,u(  •) ) ,f  (t,u,u  ( •) ) f {f  (t,u,u  (• ) ) } C (£(t,u,u(*) ) ,f  (t,u,u  ( •) ) ) iai  IQ  . 

Assume  furthermore  that  the  two  functions  f , f f Z and  are  continuous  and  bounded  on 
I x ]Rn  x Z.  Construct  the  maximal  interval  solutions  to  a,  and  a,  f to  the 
problem  (1)  , (2)  and  call  them  (v,w)  and  (v,w)  . Then  < w and  for  the  set  {u} 
of  all  solutions  u f Z of  (5) , (6)  the  inclusion 
(21)  (u)  C |v,w) 

is  true. 

Remarks : 1)  There  is  nothing  assumed  for  one  of  the  right  hand  sides  f ( If)  if 
f ^ f_, f , only  the  existence  of  that  function.  If  f for  example  is  not  continuousi  no 
solution  u of  (1),  (2)  may  exist.  Kindly  note  that  the  theorem  deals  only  with 
existing  solutions. 

2)  In  order  to  find  the  two  bounds  v,w  for  all  (in  genera]  " many)  solutions  of 
(5),  (6)  one  has  to  determine  the  maxima]  interval  solutions  of  two  coupled  systems  with 
2n  equations  each.  The  functions  v and  w arc  a "side  effect"  of  this  procedure, 
they  are  not  needed  for  the  inclusion  (21)  . They  do  have  however,  a meaning  as  "inner" 
bounds  to  (u)  in  the  sense  of  interval  mathematics. 

Proof : Let  u be  a solution  of  (5),  (6).  Then  u'  = f(t,u,u(*))  < f(t,u,u(*))  in  1^. 

Then  u < w in  I by  the  corollary  1 of  the  lemma  and  by  the  construction  of  (v.w) . 

The  inequality  u > v is  shown  similarly  which  finishes  the  proof. 
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10.  Two  examples 


1)  Let  n = 1,  a = 0, 


f(t,u,u(-))  :=  2 (/|u  (t)  | + u(T-t)) 


{f}  :=  {f 1 0 < t < 1}  . 

1 s r 

The  functional  used  is  u(T-t)  by  (9).  I do  not  know  if  the  problem  (1),  (2)  with  f 
by  (22)  can  explicitly  be  solved  for  t * 0,1. 

Since  f is  isotone  in  u(-)  the  two  inequalities  (12)  and  (13)  of  the  Lemma 


are  decoupled.  By  putting  v(t)  :=  -e,  0 < e < 1 one  gets  v(0)  = -e  < 0 = a and 

0 = v 1 (t)  < 2(/|v(t)  | + v (t  • t)  ) = 2/e  (1  - /t)  . 

Hence  by  the  Lemma  u(t)  > -e  for  any  solution  u f Z of  (1),  (2).  For  e ■+  0 one 

gets  u(t)  > 0 and  therefore  by  (22)  there  is  also  u' (t)  > 0 for  any  solution  u. 

Hence  0 < u(T-t)  < u(t)  for  any  solution,  therefore  one  can  define 

f_(t,u,u(“) ) :=  0,  f (t,  u,u  ( ■ ) ) :=  2(/|u(t)  | + u(t))  . 

— 2 

The  maximal  interval  solution  is  found  easily  as  v :=  0,  w :=  (exp  t - 1)  hence 

(23)  u(t)  e (0,  (exp  t - l)'2]  for  t > 0 , 

for  all  solutions  u of  (1),  (2)  with  (22). 

If  one  now  changes  the  functional  (9)  in  (22)  to  (10)  one  gets  the  function 


f(t,u,u(*))  :• 


2 (/|u  ( t)  | + u(t  - s))  for  0 < s < t , 


2/|u(t)  ] 


for  0 < t < s , 


where  0 < s < “>.  With  this  right  hand  side  (24)  one  gets  a whole  set  of  difference- 
differential  equations  with  retarded  argument.  The  same  ideas  as  above  give  exactly  the 
same  bound  functions  v and  w.  Hence  also  in  this  case  (23)  is  true. 

The  same  can  be  said  for  the  third  different  right  hand  side 

(25)  f(t,u,u(-))  :=  2(/|u(t)  | + (/  |u(s)  |Pds)1/P) 

0 

with  1 < p < «.  The  functional  in  this  case  is  the  Volterra  p-norm  with  the  sup  norm 
for  p » •». 
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Since  in  all  three  cases  the  bounds  v,  w are  solutions  itself  to  u’  = f,  u*  ■ 1 
one  gets  in  addition  the  optimality  condition 
(26)  v,w  e {u}  C (v,w)  . 

This  result  is  highly  surprising.  The  throe  problems  with  the  different  right  hand 
sides  (23),  (24)  and  (25)  most  certainly  have  completely  different  solutions  and  therefore 
also  different  solution  sets.  In  spite  of  this  fact  all  three  sets  have  the  same  bounds 
and  furthermore  these  bounds  are  optimal. 

This  example  shows  also,  that  it  is  very  often  simpler  to  look  for  bounds  v,w  • 
such  that  (26)  is  true  than  to  try  to  solve  the  equations. 

2)  Let  n = 2,  the  given  system  is 

{U1  = U2'  ui(0)  = “i  ' 

U2  = -V  U2(0)  = “2  • 

The  uniquely  determined  solution  is 

u1<t)  = Oj  cos  t + a2  sin  t,  u2<t)  = -a^  sin  t + a2  cos  t . 

Let  e (0,1),  e (0,1).  Then  by  the  rules  of  interval  arithmetic  (see  R.  E.  Moore 
(1966)) 

u^(t)  e (0,1)  cos  t + (0,1)  sin  t , 
u2(t)  e (-1,0)  sin  t + (0,1)  cos  t . 

The  set  of  solutions  is  hatched  in  Figure  1.  In  the  picture  also  the  "main"  solution 
for  a\  = a2  = *S  shown*  T^e  extended  system  (19)  reads  here  as 


vi  - 

V 

v2  (0) 

= o , 

"a" 

"Wl' 

v2  (0) 

» 0 , 

V3  " 

W2' 

v3(0) 

- 1 . 

V4  “ 

-Vl' 

v (0) 

4 

= 1 . 

The  solution  to  this  system  is  unique,  hence  the  maximal  interval  solution 


$ 


1 


2v^(t)  = -efc  + sin  t + cos  t , 

2v2(t)  = -e1  + cos  t - sin  t , 

2w^ (t)  = efc  + sin  t + cos  t , 

2v<2  (t)  = et  + cos  t - sin  t . 

Hence  the  functions  v,  w "back  away"  from  the  "main"  solution  (see  Figure  1)  as  fast 
t 

as  e /2  to  below  and  to  above. 

For  t = 2n  one  sees  {u,  (2u)  } , {u_  (2ir)  } C [0,1].  But  w,  (2fl)  -v.(2ti)  = e2"  ='53&.4  ... 

1 2 — k — k 

for  k = 1,2.  The  real  set  of  solutions  {u}  is  therefore  surpassed  at  t = 2v  by 
[v,w]  by  a factor  of  more  than  500  and  this  grows  rapidly  worse  for  larger  values  of  t ! ! ! 
This  example  was  first  discussed  by  R.  E.  Moore  (1966)  by  using  geometrical  reasoning. 


1 


11 . Optimal  bounds  for  the  set  of  solutions;  maximal  and  minimal  solutions 

By  section  9 one  sees  that  it  is  enough  to  restrict  the  survey  to  the  solution  of 
two  real  functional  differential  equations  if  one  wishes  to  bound  the  solutions  of  sets 
of  such  equations.  Hence  the  following  optimality  considerations  are  given  only  for 
systems  of  the  type  (1) . Since  it  is  no  aggravation  the  initial  inclusion  (6)  has  however, 
been  used  instead  of  the  initial  condition  (2).  It  is  always  assumed  that  (15)  is  true. 

As  is  shown  by  the  second  example  of  section  10  the  maximal  interval  solution  (v,w) 
gives  in  general  not  optimal  bounds  to  the  set  of  solutions  (u).  In  this  and  in  the 
next  section  classes  of  problems  will  be  given  such  that  there  is  optimality  either  in 
the  sense  of  (26)  or  at  least  of  (8) . 

Let  f e M.  Assume  furthermore  that  all  functions  f,  (t,u,  ,.u,v(*))  are  uncondi- 

k k k 

tionally  partially  isotone  with  respect  to  any  of  the  components  of  aTxl  v(-). 

This  is  called  "quasimonotone  increasing"  by  W.  Walter  (1970)  in  the  case  of  differential 
equations.  With  this  condition  the  equations  for  v and  w in  (16)  and  (19)  are 
decoupled  from  each  other.  Therefore  the  functions  v and  w are  even  solutions  of 
(1) . In  the  case  of  differential  equations  these  are  the  well  known  minimal  and  maximal 


12.  Further  cases  with  optimal  bounds 

If  one  changes  in  the  system  (26)  the  sign  in  the  second  equation  one  gets 

C U*  = 2/fuj".  11^(0)  = 0 , 


1 u*  = -2/|u1|,  u2(0)  = 0 . 

Now  the  right  hand  sides  are  not  anymore  "quasimonotone  increasing".  The  set  of  all 

a 2*  2 

solutions  can  be  described  quite  easily,  one  finds  { u^ ( t ) } C (0,t  ),  {u2<t)}  C (0,t  ). 
The  maximal  interval  solution  of  (19)  , (20)  produces  the  bounds 

Vj(t)  :=  0,  v2(t)  :=  -t2,  w^(t)  :=  t2,  w2  :=  0 . 

They  are  again  optimal  bounds.  Opposite  to  the  results  of  section  11  the  functions  v 

and  w are  not  anymore  solutions  of  (1)  (but  certainly  of  (19)). 

By  inspection  one  sees  however,  that  the  "crossed"  couples  (v  w2>  and  (w^,v2) 

each  are  a solution  to  (1).  This  is  responsible  for  the  optimality.  Because  if  v 

and  w are  at  least  componentwise  solutions  of  (1)  then  there  exist  no  smaller  intervals 
[v,wj  C fv,w]  such  that  (4)  is  valid.  Hence  in  this  case  (8)  is  true  which  means 
optimality. 

The  classes  of  problems  (1),  (6)  considered  in  this  section  are  extensions  of  this 
example . 

Definition  (monotonicity  matrices) : Let  f e M. 

Define  for  i = l(l)n 

3ii  !=  1 ' 


r° 

if 

f . 

1 

does  not  depend  upon 

Uk  ' 

aiK  < +1 

if 

fi 

depends  isotone  upon 

\ ' 

u 

if 

£i 

depends  anti  tone  upon 

Uk 

for 

i * k - l(l)n  , 

0 

if 

f . 

1 

does  not  depend  upon  u^ ( ■ ) , 

+1 

if 

fi 

depends  isotone  upon  u^ ( • ) , 

-1 

if 

f . 

1 

depends  antitone  upon 

for 

k = 1 (l)n  . 
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The  matrices  A - (a^)  a**i  B - (b^)  are  called  the  monotonic  matrices  to  f. 
Definition  (monotonicity  condition  (M) ) : Let  f t M. 

* i 

Assume  the  existence  of  associate  matrices  A'  = (a'  ),  B'  = (b’  ) to  the 

ik  lk 


monotonicity  matrices 

A and  B 

such 

that 

for  all 

i,k  = 1 (1) n 

f 

aik 

e { + 1, 

-1). 

b!  * 
ik 

{+1,-1}  , 

aik 

= aik 

for 

a.,  * 0 
ik 

! 

(28) 

< 

blk 

= bik 

for 

b.,  * 0 
ik 

$ 

aIk 

= bik 

$ 

k. 

aik 

= aii 

• aik 

• 

Then  f is  said  to  satisfy  the  condition  (M) . 
Remark : By  (28)  one  sees  immediately 


ik 


= a't  • a'^  for  all  l = l(l)n  . 


Ik 


Theorem: 


Assume  (15).  Let  f t M be  continuous  and  bounded  on  X x jRn  x Z.  Assume  that 
f satisfies  the  condition  (M) . Let  lv,w]  be  the  maximal  interval  solution  to  (1) , (6)  . 

Then  (v,wj  is  even  the  interval  hull  of  the  set  of  all  solutions  {u}  of  (1),  (6), 
i.e.  (8)  is  true. 

Proof : By  construction  of  [v,w]  the  inclusion  (18)  is  true.  Define  the  two  function 

vectors  p = (p  ,p P > = p(t),  q = (q  ) = q(t)  by 

i ^ n i / n 


- k 

for 

aik 

* +1  , 

"k 

for 

“lk 

= -1  , 

Wk 

for 

aik 

= +1  , 

2^k 

for 

aik 

= -1  . 

Certainly  p^(0)  * (a)  and  (0)  t (a).  Moreover  the  vectors  p and  q are  by 
construction  and  by  (28)  both  solutions  of  (1).  Hence  v and  w are  componentwise 
composed  of  solutions  of  the  problem  (1),  (20).  Hence  (8)  is  true. 


•) 


These  need  not  be  uniquely  determined. 
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